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Abstract 

m . 

I> ' We consider a preferred- frame bimetric theory in which the scalar gravitational 

field both influences the metric and has direct dynamical effects. A modified 
^ '. version ("v2") is built, by assuming now a locally-isotropic dilation of physically 

CN ' measured distances, as compared with distances evaluated with the Euclidean 

^ ■ space metric. The dynamical equations stay unchanged: they are based on a 

^3 ■ consistent formulation of Newton's second law in a curved space-time. To obtain 

I a local conservation equation for energy with the new metric, the equation for the 

O^' scalar field is modified: now its l.h.s. is the flat wave operator. Fluid dynamics 

5-H . is formulated and the asymptotic scheme of post-Newtonian approximation is 

adapted to v2. The latter also explains the gravitational effects on light rays, as 
' did the former version (vl). The violation of the weak equivalence principle found 

^ . for gravitationally-active bodies at the point-particle limit, which discarded vl, 

^ I is proved to not exist in v2. Thus that violation was indeed due to the anisotropy 

■ ■ ' of the space metric assumed in vl. 
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1 Introduction 

The universality of gravitation, i.e., the fact that "all bodies fall in the same way in a 
gravitational field," is a distinctive feature of the gravity interaction. It is also known 
as the "weak equivalence principle" (WEP), the equivalence being between the gravita- 
tional and inertial effects: indeed, the fact that the gravitational acceleration of a small 
body is the same independently of its composition and its mass, allows to incorporate 
in it the acceleration due to the motion of the arbitrary reference frame through an 
inertial frame. Following Galileo, Newton, and Eotvos, the empirical validity of the 
WEP has been checked with increasing precisions jTj. Therefore, the WEP is involved 
in the very construction of relativistic theories of gravitation, in particular in the con- 
struction of Einstein's general relativity (GR), which takes benefit of the universality 
of gravity to geometrize it. More precisely, in the construction of a relativistic theory. 
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one ensures the validity of the WEP for test particles, which, by definition, do not 
influence the gravitational field. The geometrization operated by GR is one way to do 
this. Another way is to formulate the dynamics by an extension of Newton's second 
law, the gravitational force being mg with m the (velocity-dependent) inertial mass 
and g a (theory-dependent) gravity acceleration 2J. We would like to emphasize that 
Einstein's geometrized dynamics, in which test particles follow geodesic lines of the 
space-time metric, can be rewritten in that way [2]. 

A theory of gravitation has been built, in which one starts from a simple form for 
the field g, suggested by a heuristic interpretation of gravity as Archimedes' thrust 
due to the space variation of the "ether pressure" pe, and in which pe also determines 
the space-time metric. One obtains a scalar bimetric theory with a preferred reference 
frame. This theory is summarized in Ref. |3, or in more detail in Ref. It admits 
a "physically observable preferred foliation" [S], in short a detectable ether. Yet, apart 
from the WEP violation to be discussed below, it seems to agree with observations 
E]. In particular, the ether of that theory could be indeed detected by adjusting 
on astrodynamical observations the equations of celestial mechanics that are valid in 
the theory One of the motivations to investigate a preferred-frame theory is that 
the existence of a preferred reference frame, involving that of a preferred time, may 
be regarded as a possible way to make quantum theory and gravitation theory match 
together [Hj. It is interesting to note that, currently, a detectable ether is advocated 
already in the absence of gravitation by Selleri, in the framework of a theory close to 
special relativity, but based on an absolute simultaneity 0. 

As we mentioned, the WEP is valid for any test particle, by construction — in the 
investigated scalar theory just as well as in GR. But any real body, however small it 
may be, must influence the gravitational field, so that one can a priori expect the oc- 
currence of self-accelerations. The latter ones are excluded from Newton's theory by 
the actio-reactio principle, but this principle cannot even be formulated in a nonlinear 
theory — as are most relativistic theories. Moreover, the mass-energy equivalence, which 
has to be accounted for by a relativistic theory, implies that the rest-mass, kinetic, and 
gravitational energies of the test body all influence the gravitational field. This means 
that the internal structure of the test body may a priori be expected to influence its 
motion, which would be a violation of the WEP. Hence, the validity of the WEP for test 
particles is far from ensuring that this principle applies to real bodies. In GR, studies of 
the self-force acting on an extended body can be found, e.g. 013121110], but they have 
been based on simplifying assumptions such as that of a black-hole body [ZllHllini, or on 
hnearized GR 0. Such assumptions do not seem very appropriate to check whether a 
violation of the WEP might be predicted by GR in the real world, say for an asteroid or 
a spacecraft in the solar system. It is well-known, since the work of Nordtvedt fn\ IT^ , 
that the WEP may be violated for extended bodies of a finite mass in some relativistic 
theories of gravity (see also Will El)- However, the kind of violation of the WEP 
which we have in mind is a more severe one, that occurs even in the limit in which the 
size of the extended body shrinks to zero P^ . 
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In the above-mentioned scalar theory ^ , it has been proposed a rigorous framework 
for the study of weakly-gravitating systems ^H] , in conformity with the asymptotic ap- 
proximation schemes which are currently used in applied mathematics, and this scheme 
has been developed up to the equations of motion of the mass centers of perfect-fluid 
bodies [T71 118j. This "asymptotic" scheme of post-Newtonian approximation (PNA) is 
based on a one-parameter family of initial conditions, defining a family of gravitating 
systems. ^ It is hence different from the "standard" PNA proposed for GR by Fock 
PU] and Chandrasekhar which is at the basis of a significant part of the subse- 
quent work on relativistic celestial mechanics in GR, and in which no such family is 
considered. The asymptotic PNA predicts that the internal structure of the bodies def- 
initely influences their motion in a weakly-gravitating system such as our solar system, 
at least for the scalar theory ^H]. Moreover, by considering a family of PN systems 
which are identical but for the size of one of the bodies, which is a small parameter ^, 
it has been possible to make a rigorous study of the point-particle limit. It has thus 
been found jTH] that the internal structure of a body does influence its post-Newtonian 
acceleration even at the point-particle limit. It has also been investigated the particular 
case where, apart from the small body, there is just one massive body, which is static 
and spherically symmetric (SSS). In that case, the PN equation of motion of the mass 
center of the small body is identical to the PN equation of motion of a test particle 
in the corresponding SSS field, apart from one structure- dependent term. These results 
show a patent violation of the WEP, and one whose magnitude is likely to discard the 
initial version of the theory ^21- The specific reason which makes the WEP violation 
actually occur in the point-particle limit of the PN equations of motion has been iden- 
tified: it is the anisotropy of the space metric (in the preferred reference frame of the 
theory). More precisely, it was assumed that there is a gravitational contraction of 
physical objects, only in the direction of the gravity acceleration g (see Ref. j22] and 
references therein), thus making the PN spatial metric depend on the derivatives U^i of 
the Newtonian potential U. The structure-dependence (hence the non-uniqueness) of 
the point-particle limit comes then [15J from the fact that the self part of the second 
derivatives U^ij is order zero in ^. ^ 

Therefore, it has been began |22j to investigate the case where one assumes a locally 
isotropic gravitational contraction, which case is nearly as natural as the anisotropic 

^ A scheme similar to our "asymptotic" scheme has been previously proposed for GR by Futamase 
& Schutz PH], though it is based on a particular initial condition as to the space metric and its time 
derivative, which enforces the spatial isotropy of the metric. In any case, the work |19) is restricted to 
the local equations, and this in a form which is not very explicit. 

^ The PN spatial metric valid for the standard form of Schwarzschild's metric of GR also depends 
on the derivatives U,i (with specifically U = GM/r in the case of Schwarzschild's metric). Hence, the 
dependence of the spatial metric on the U^i 's should be found for general situations with "Schwarzschild- 
like" gauges, i.e., with gauge conditions under which the standard form of Schwarzschild's metric is 
the unique SSS solution of the Einstein equations with Newtonian behaviour at spatial infinity. (It is 
not difficult to exhibit such gauges.) One may then ask ^2] whether a similar violation of the WEP 
might appear also in GR, in such Schwarzschild-like gauges. This would be more difficult to check, 
due to the greater complexity of GR as compared with the present scalar theory, in particular due to 
the necessity of satisfying constraint equations. 
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case according to the heuristic concept that led to the scalar theory, i.e., gravity seen 
as due to the heterogeneity of the field of "ether pressure" pe- It has been shown that 
there is some freedom left for the equation governing the field pe- The aim of the 
present paper is to propose a definite equation for the scalar gravitational field, leading 
to an exact local conservation law for the total (material plus gravitational) energy; 
to investigate the main features of this new version of the theory; and to show that it 
does eliminate the WEP violation which has been found with the former version. The 
next Section summarizes the basic concepts of the theory, which remain true for the 
new version to be built in this paper. Section |H1 precises the equation for the scalar 
field, in connection with the necessity of obtaining an exact energy conservation. The 
equations of motion of a perfect fiuid are obtained in Sect. IH Sectional is devoted to 
the post-Newtonian approximation and shows, in particular, that this theory predicts 
the same effects on light rays as the standard effects known in GR. The PN equations 
of motion of the mass centers are derived in Sect. IHl The point-particle limit is taken 
in Sect. [71 and our conclusion makes Sect. |H1 

2 Basic concepts 

The idea according to which gravity would be Archimedes' thrust due to the macroscopic 
part of the pressure gradient in a fiuid "ether," and the necessity to account for special 
relativity, lead to set a few basic assumptions [22] , which we state directly here. 

2.1 The preferred reference frame and the metric 

The space-time manifold V is assumed to be the product R x M, where M is the preferred 
reference body, which plays the role of Newton's absolute space. The equations of the 
theory are primarily written in the preferred reference frame E associated with that 
body. ^ Specifically, it is assumed that the preferred reference body M is endowed 
with a time- invariant Euclidean metric g^, with respect to which M is thus a rigid 
body, which is assumed to fill the Euclidean space. On the "time" component of the 

■^A reference frame is for us essentially a (three-dimensional) reference body, plus a notion of time. 
Let us begin with the viewpoint of "space plus time," which is admissible once we take the space-time 
to be a product, V = R x M — the "time" T of an event X = (T, x) S V being thus the canonical 
projection of X into R. From this viewpoint, a general reference frame can be defined as a time- 
dependent diffeomorphism ipr of the space M onto itself. Consider the trajectories T i— -0t(x), each 
for a fixed x S M. As x varies in M, the set of these trajectories defines a deformable body N, which 
is uniquely associated with the reference frame (V'r)TGR- Thus, x' = ^/'t(x) describes the motion of 
N relative to M. For a general ipr, this body N will indeed be deformed. The most obvious reference 
frame is yet that one which is associated with M itself, thus: VT, ipT = IdM- This is our preferred 
frame denoted by E. From the viewpoint of "space-time," which is more general and which is hence 
also admissible in an ether theory, a reference frame is defined by a three-dimensional congruence of 
world lines, which defines a "body" (as we say) or "reference fluid" [23]. Thus the trajectories of the 
"space plus time" description are replaced by world lines: s i-^ X{s). Among systems of space-time 
coordinates: X i-^ X^{^) — {i^ — 0: 3), coordinate systems adaptedto a given reference frame are 
such that each world line X{s) belonging to the corresponding body has constant space coordinates: 
Vs, x'iXis)) ^y' = Constant {i = 1, 2, 3). 
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space-time, we have the one-dimensional Euchdean metric. Combining these metrics 
on the component spaces, we get a Lorentzian metric: the square scalar product of an 
arbitrary 4-vector U = (f/°,u), with u an arbitrary spatial vector {i.e., formally, an 
element of the tangent space TMx to M at some x G M), is 

7°(U,U) = (t/°r-g°(u,u). (1) 

This is the "background metric," which should determine the proper time along a 
trajectory, if it were not for the metrical effects of gravity. For this to be true, it is 
necessary that the canonical projection of an event X G V gives in fact = cT, rather 
than T, where T is indeed a time (called the "absolute time") and c is a constant — the 
velocity of light, as measured with "physical" space and time standards. If on M we 
take Cartesian coordinates (x*) for g^, then the space-time coordinates (x^) = (x*)) 
are Galilean coordinates for 7° [i.e., in that coordinates, (7°j,) = diag(l, —1, —1, — 1)], 
which are adapted to the preferred frame E. The gravitational field is a scalar field 
Pe, the "ether pressure," which determines the field of "ether density" pe from the 
barotropic relationship = Pe{Pe)- Gravity has metrical effects which occur through 
the ratio 

/?(T,x)^Pe(T,x)/pr(T)<l, (2) 

with 

pT (T) = Sup^gMPe (T, x) . (3) 

More precisely, the "physical" space-time metric 7, that one which more directly ex- 
presses space and time measurements, is related to the background metric 7° by a dila- 
tion of time standards and a contraction of physical objects, both in the ratio j3. Thus, 
in any coordinates (y^) which are adapted to the frame E and such that = = cT, 
the line element of 7 writes 

ds' = ^.Ay^dy'' = /32(dy°)2 - g^,dy'dy^, (4) 

where g is the physical space metric in the frame E. This equation implies that the 
scalar (3 can be more operationally defined as 

P^{loof" (5) 

(in any coordinates {y^) adapted to the frame E, and such that y^ = cT). Metric g is 
related to the Euclidean metric g^ by the assumed gravitational contraction of objects, 
including space standards — hence the dilation of measured distances. This effect was 
formerly assumed to occur only in the direction of gradpe (see e.g. Ref. jl]), but it 
is now assumed isotropic. This means that the following relation is now assumed [221 
between the spatial metrics and g: 

9 = rV- (6) 
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2.2 Dynamical equations 

In addition to its metrical effects, the gravitational field Pe has also direct dynamical 
effects, namely it produces a gravity acceleration 

g ^ (7) 

Pe 

[The index g means that the gradient operator refers to the physical, Riemannian metric 
g, i.e., (gradg0)(x) is a spatial vector (an element of TMx) with components (gradg0)* = 
g^-'(pj.] We shall soon assume that the Pe-Pe relationship is simply p^ = c^Pe- In that 
case, we get from Q: 

grad /? 

g = -c'— (8) 

Equation ((7j) is a fundamental one according to the concept of the theory. However, to 
use this "gravity acceleration," we must define the dynamics by an extension of Newton's 
second law, as announced in the Introduction. This has been done in detail in Refs. 
[211211123 (see Ref. |1] for a summary). The dynamical equations stay unchanged in 
the new version of the theory. Therefore, only a synopsis of that alternative dynamics 
will be given here. Dynamics of a test particle is governed by a "relativistic" extension 
of Newton's second law: 

. . DP 

Fo + m{v)g = (9) 



''X 



where Fq is the non-gravitational {e.g. electromagnetic) force, v = d'^/dty^ the velocity 
and V = gf(v,v)^/^ its modulus, both being measured by physical clocks and rods 
of observers bound to the preferred frame E; m{v) = m(0)7„ is the relativistic mass 
(7„ = (1— f ^/c^)^^/^ is the Lorentz factor); P = m(f )v is the momentum; tx is the "local 
time" in E, measured by a clock at the fixed point x in the frame E, that momentarily 
coincides with the position of the test particle: from (Q, we have 

dt^/dT = p{T,x); (10) 

and Dw/D^ is the appropriate "time-" derivative of a vector w(^) in a manifold M 
endowed with a "time-" dependent metric (21 121! • static case (/?,o = 0) with 

Fq = 0, that dynamics implies Einstein's motion along geodesies of the curved space- 
time metric 7 (see Ref. 0). ^ 

Dynamics is also defined for continuous media. For a dust, we may apply Q point- 
wise and this implies PHI the following equation: 

T;;. = b„ (11) 

^ /3 being defined by Eq. jSJl, postulating Eq. ^ for vector g is equivalent, under natural require- 
ments, to asking geodesic motion in the static case 0. On the other hand, if one adds to the r.h.s. of 
Eq. ^ a complementary term involving the velocity of the test particle and the time-variation of the 
space metric, then Eq. implies geodesic motion in the general case, thus one deduces Einstein's 
geodesic motion from a relativistic extension of Newton's second law |2j. 
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where T is the energy-momentum tensor of matter and nongravitational fields, and 
where is defined by 



(Indices are raised and lowered with metric 7, unless mentioned otherwise. Semicolon 
means covariant differentiation using the Christoffel connection associated with metric 
7. Note that in GR, in contrast, we have 6^ = in The universality of gravity 

means that Eq. (fTT|) with definition (fT^ must remain true for any material medium. 
Thus, the dynamics of a test particle, as well as the dynamical field equation for a 
continuous medium, exactly obey the WEP, just as it is the case in GR. Equations 
flTT|) - (fT^ are valid in any coordinates {y'^) which are adapted to the frame E and such 
that y^ = 0(T) for some function (p. 

3 Scalar field equation and energy conservation 

3.1 Semi- heuristic constraints on the scalar field equation 

Equation ((Zj) for the gravity acceleration expresses the idea according to which gravity 
would be Archimedes' thrust in Romani's "constitutive ether" [2E], i-G. a space-filling 
perfect fiuid, of which any matter particle should be a mere local organization. The 
equation for the field Pe in ((Tj) should fulfil the following conditions [221]: 

i) Newton's gravity, because it propagates instantaneously, should correspond to the 
limiting case of an incompressible ether (pe = Constant). 

ii) In the real case, the ether should have a compressibility K = 1/c^, so that the 
velocity of the pressure waves, Ce = (dpe/dpe)^^^ (beyond which velocity the material 
particles, seen as fiows in the universal fiuid, should be destroyed by shock waves), 
coincide with the relativistic upper limit c. (This constraint leads to Pe = c^Pe, as 
assumed in advance for Eq. (jH)).) Thus, in the real case, the equation for pe should 
be a kind of nonlinear wave equation, the nonlinearity arising from the fact that the 
physical space-time metric 7 is determined by the field pe itself. It is indeed the physical 
metric, not the background metric, which is directly relevant here, because the rela- 
tivistic upper limit c applies to velocities measured with the local, physical instruments. 

We note that, from Eqs. (jH) and (jH)), the physical metric 7 is nearly equal to the 
given Galilean metric 7° if and only if /? ^ 1, which, by the definition means a 
quasi-incompressible fiow, pe ~ Constant if moreover we consider a "non-cosmological 
time scale" so that pf ~ Constant. (Due to Eq. (jSJ, it also implies that the gravita- 
tional field is weak, i.e. the gravity acceleration is "small.") Thus, in the limit /5 — > 1 
with ~ Constant, the nonlinearity of the wave equation must evanesce while si- 
multaneously the ether compressibility becomes very small and the metric becomes 
close to Galilean. Now Newton's gravity does correspond to a Galilean metric and is 
characterized by Poisson's equation for the field g: 




(12) 



divgo g = —AnGp, 



(13) 
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with G the gravitational constant and p the Newtonian mass density. Hence, for an 
incompressible ether (pe = Constant), in which the gravity acceleration is defined by 
Eq. (|7j) with g = g^, Newton's gravity is exactly equivalent to the following equation 
for the "ether pressure" Pe- 

Agope = 4:7iGppe. (14) 

Therefore, the conditions i) and ii) suggest to admit that the equation for the field 
Pe has the form j22] 

AgPe + (time derivatives) = AnGapeF^P), (15) 



F{P) ^ 1 as /3 ^ 1, (16) 

the "time derivatives" term being such that, in the appropriate ( "post-Minkowskian" ) 
limit 15, involving the condition (3^1, the operator on the l.h.s. becomes equivalent 
to the usual (d'Alembert) wave operator. On the r.h.s. of (|TH|l . a denotes the relevant 
mass-energy density, which shall have to be defined in terms of the energy-momentum 
tensor T (that we take in mass units), and which will reduce to the Newtonian (rest- 
mass) density p in the nonrelativistic limit. It is worth noting that the equation assumed 
^ for the previous (anisotropic-metric) version of the theory is indeed a special case of 
Eq. (fT^ Let us thus start from (fT^ and (fTB|). or equivalently from (fTH|) and 

AgPe + (time derivatives) = — —apeF{/3). (17) 

3.2 The equation for the scalar gravitational field and the en- 
ergy equation 

As we mentioned, the physical metric is directly relevant to translate the heuristic con- 
siderations on the "gravitational ether" into restrictions imposed on the field equation. 
However, the latter would obviously be more tractable if it could be rewritten in terms 
of the flat background metric. Let us try to impose this condition. To this aim, we 
evaluate the spatial term Ag p^ on the l.h.s. of (fT7|) . From the definition 

Ag<p = divg(gradg0) = {V9 9''<Pj),, {9 = det{g,,), (g^^) = (g,,)-') , (18) 

and since, by (jHI), we have 

g = p-'g' (^°^det(4)), (19) 
it follows that, in Cartesian coordinates for g^, it holds 

^0Pe = P'(¥) ■ (20) 



Accounting for the definition of /3 from pe (Eq. (0)), we get: 

AgPe = pTf3' (^^) = pTf3' Ago (Log/?) . (21) 
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Hence, if the sought-for finalization of Eq. ()17|) is to be reexpressed nicely in terms of 
the flat metric, the field variable should be 



i/j = -Log/3. (22) 
(The minus sign is chosen so that ip > 0, see Eq. 1^.) 

On the other hand, we impose on the field equation the additional condition that 
an exact local conservation law must be found with some consistent definition of the 
energy. The latter should be the sum of a material energy and a gravitational energy. 
The energy conservation law should be obtained by rewriting the time component of 
the dynamical equation, (fTT| with the definition ()12j) . as a zero 4-divergence. For any 
metric having the form (^, and independently of any restriction on the equation for 
the scalar field, the following identity: 

TH;. = ^ (v^^;),. - (7 = det(7,.)) (23) 

allows us to rewrite the time component of as 

(v^^c?), + (v^^o),o = v^moT°° = a. (24) 

Thus, a must be a 4-divergence by virtue of the equation for the scalar field. Using 
now the specific form (jH)) assumed for the space metric, we have Eq. (fT^ and get from 
0: 

-7 = A = r'^7°, (25) 
so that, in Cartesian coordinates for and with = cT, 

a=(Log/3)oTOO^-^,oTOO. (26) 



Together with (j2Tj) . Eq. (j26|) makes it obvious that the relevant scalar field is indeed 
■ip = —Log/3. We note the identities 

V'.o = (V-.o - \ (V^jV^,i),o ' V-.o V'.o.o = \ (^,o) ,0 • (27) 

From these, and from ()26j] . it follows that, if we take the gravitational source o to be 
the energy component T°°, and if we simply postulate for -0 the flat wave equation: 

U^ = i)^^^^- ^g.i) = ^a {a = T'^\ x'' = cT), (28) 



c 

then we indeed obtain a as a 4-divergence: 

« = ^ { -|5o [^'o + (gradV^)'] + div (V^,ogradV^) } . (29) 

(Henceforth, div, grad and also A shall be the standard operators defined with the 
Euclidean metric g^.) More complicated time derivatives in the field equation could 
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also provide a conservation equation, but the spatial term is more or less imposed to 
be /S.ip by (PT|) . while the source a should be due to (pUj). We have currently few 
constraints on the time-derivative part of the equation for the scalar gravitational field. 
Hence, Occam 's razor leads us to state as the equation for the scalar field. (This 
corresponds to F{P) = in Eq. (fT7|) .) 

Thus, by assuming the validity of Eq. (j^ . we rewrite a [Eq. (j^ ] as (j^ . Hence, 
becomes the following local conservation equation for the energy: ^ 

9o(£m + eg) + div($^ + $g) = 0, (30) 

where the material and gravitational energy densities are given (in mass units) by: 

^ T°°, = ^ W + (grad^)'] , (31) 

and the corresponding fluxes are: 

$„ = (T°^), $g ^ --^ (^^ograd^/.) . (32) 

The scalar field equation ()28j) and the conservation equation fl3()j] are valid in any 
coordinates {y'^) adapted to the preferred frame and such that = x° = cT. We 
note in particular that, although the d'Alembert operator □ is generally-covariant, Eq. 
(I2H1) does not admit a change in the time coordinate ?/°, because ip = ~Log^/7oo and 
cr = behave differently under a change y'^ = 0(?/°). ^ Thus, Eq. (j^Hj) admits only 
purely spatial coordinate changes, consistently with the preferred-frame character of 
the theory. (It is recalled at § 15.31 how to cope with this character, on the example of 
the effects on light rays; see at the end of Ref. [3 for the case with celestial mechanics.) 



3.3 Comments on the balance equation for the spatial momen- 
tum 

Similarly, let us rewrite the spatial component of the equation of motion of a continuum 
fjllj) in terms of the scalar field (j^^ . using the explicit form (0])-® of the metric, and 
the identity 

T^'^ = iV^T'^") + K^.T-^' (33) 

(where the T'j^'s are the Christoffel symbols of metric 7). Adopting Galilean coordinates 
(x^) for the flat metric 7° henceforth, we find after an easy algebra: 

(g4V- jij^ ^ ^ (^g4V. yio^ ^ _ g4V. + ^^^jio^j ^ (34) 

5 We use the fact that, from Q and jSU, yf^Tl^ = P^'^T^ = T"" and V^^o = (^'"^T^o = 
in Cartesian coordinates, so that H3U|) with (|31|l and (|32|l apply then — but these are space-covariant 
equations. 

^ However, we note also that a mere change in the time unit, T' — aT, does not affect the time 
coordinate a;" = cT (since c becomes c' ~ c/a), hence leaves Eqs. (|28|) and (|3U|I invariant. 
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An identity similar to (P7|) (with ip^i instead of ifj^) allows to get the r.h.s. as a 4- 
divergence, using the scalar field equation (j2HI)- Due to the remaining source term on 
the l.h.s., it is in general not possible to rewrite ()34|) as a zero fiat 4-divergence. I.e., 
there is no local conservation equation for the total (material plus gravitational) mo- 
mentum in this theory. 

In contrast, in Lagrangian-based relativistic theories of gravitation, e.g. in GR, there 
is a local conservation law (or something that looks like that) for the total momentum, 
which is the sum of the local momentum of matter and the local (pseudo-)momentum 
of the gravitational field. (The meaning of the latter decomposition and of its coor- 
dinate dependence is clearer in the teleparallel equivalent of GR t2Zi-) In some cases, 
characterized by a sufficient fall-off of at spatial infinity, the global value (i.e. the 
space integral) of the total momentum is then conserved. However, this does not mean 
that the global value of the momentum of matter is then conserved: in fact, it precludes 
this, unless the global momentum of the gravitational field is separately constant — but 
this occurs only when the gravitational field is constant, thus when there is no motion of 
matter. Therefore, the situation is not so much different in the investigated theory and 
in Lagrangian relativistic theories: in both kinds of theories, the global momentum of 
matter is in general not conserved, unless there is just one body in equilibrium — the lat- 
ter case is, of course, possible also in the investigated theory. [Assume time-independent 
fields in Eqs. fl28j) . (jHn|l and (jSl}-] In Newton's theory, in which there is a local conser- 
vation for the total momentum, the global momentum of matter is conserved, however. 
This is because the global value of the gravitational momentum turns out to be zero 
in Newton's theory j21j. (The physical reason for this is that there is no gravitational 
radiation in Newton's theory.) The fact that, in contrast, the momentum of matter is 
in general not conserved in relativistic theories of gravitation, is related to the generic 
presence of self-accelerations (or self- forces) in these theories (including GR), already 
mentioned in the Introduction. 

4 Equations of motion and matter production for a 
perfect fluid 

In most applications of a "relativistic" theory of gravitation, it is enough (and it is 
indeed usual [201 1211 I2S1 ISOl El) to consider a perfect fiuid, because: i) the stress 
tensor, i.e., the spatial part of tensor T, has normally a non-spherical part small enough 
that the latter does not bring significant post-Newtonian (PN) corrections; and ii) the 
motion of astronomical bodies can be described as approximately rigid (here also, it 
is an even better approximation if one assumes this only at the stage of calculating 
the PN corrections), in which case a viscosity has no effect. For a perfect fluid, with 
its well-known expression for tensor T PO], depending on the pressure p, the proper 

^ In fact, it does not seem completely clear what should be the physically motivated conditions 
ensuring the sufficient decrease at infinity for 0, due to the fact that one has to account for gravitational 
radiation: see e.g. Stephani 
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density of rest mass p*, the density of elastic energy per unit rest mass 11, and the 
velocity u = dx/dT = /?v, we introduce the field variable 



(35) 



/32 c2 /32 

and rewrite the equation of motion ()34|1 and the energy conservation (jHIH) respectively 

{Ou') y + {eu'u^) . - tfj^T Ou' - tfj^i Ou^u^ = c^%l)^i a + e^'^ (p t/;,^ - p,^) (36) 



as: 



and 

[e-'^^e) ^ + [e-'^eu^) ^ = -tT a+\ (e'^p) ^ . (37) 



As it has been discussed in detail in Ref. [25j, the exact energy conservation of the 
scalar theory precludes in general an exact conservation of (rest-) mass. There, it has 
been shown that, already for a perfect and isentropic fluid, the general form of the 
equations of motion for a continuum implies a reversible creation/ destruction of matter 
in a variable gravitational field. Let us note U the 4- velocity, with f/^ = dy'^/ds. One 
finds that the rate of creation/destruction is jSSj: 

(38) 

This equation holds true independently of the scalar field equation and the specific form 
of the space metric [221 • With the new form (jHI) assumed for the space metric, we get: 

$ = -6^, (39) 

which is three times the rate found with the formerly-assumed anisotropic space metric 
Of course, the actual values of P and jS^ in a given physical situation may depend 
on the theory. However, anticipating over the next Section, we can write down the 
weak-field approximations of (3 and its relative rate as: 

with U the Newtonian potential, whose time-derivative has to be taken in the preferred 
frame. This, indeed, gives three times the former weak-field prediction for p, namely it 
gives 

P-^Of^ (41) 

which remains extremely small in usual conditions (but would be significant inside 
stars): the relative creation rate p/p would be ca. 10~^^s~^ at or near the surface of 
the Earth, if the "absolute" velocity of the Earth is taken to be 300km.s~^. {See Ref. 
Sect. 4.3. Note that equal amounts of mass would be destroyed and created at 
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opposite positions on the Earth, Eq. (4.22) there.} Mass conservation is far to have 
been checked to this accuracy. Note that matter creation is being actively investigated 
in cosmological hterature, see e.g. jHH ISS EH ESI EE] • If niass non-conservation is to 
occur in "cosmological" conditions, it must exist in nature, and so possibly (in minute, 
not yet observable quantities) in today's world. 

5 Post-Newtonian approximation (PNA) 

5.1 Definition of the asymptotic scheme of PNA 

The purpose of the post-Newtonian approximation is to obtain asymptotic expansions of 
the fields as functions of a relevant field-strength parameter A, and to deduce expanded 
equations (which are much more tractable than the original equations) by inserting the 
expansions into the field equations. To do this in a mathematically meaningful way, 
it is necessary that one can make A tend towards zero, hence one must (conceptually) 
associate to the given gravitating system S a family (Sa) of systems, i.e., a family of 
solution fields of the system of equations. This family has to be defined by a family 
of boundary conditions — initial conditions for that matter, because here gravitation 
propagates with a finite velocity. The system of interest, S, must itself correspond to 
a small value Aq of the parameter, thus "justifying" to use the asymptotic expansions 
for that value. The definition of the family involves two conditions which make this 
family represent the Newtonian limit: as A ^ 0, i) the physical metric 7*^^^ must tend 
towards the flat metric 7", and ii) the fields must become equivalent to "corresponding" 
Newtonian fields. The first condition is easy to be explicited in a scalar theory, in which 
the relation between 7 and 7*^ depends only on the scalar gravitational field. Condition 
ii) asks for two preliminaries: a) that one disposes of a relevant family of Newtonian 
systems, for comparison, and b) that one is able to define a natural equivalent of the 
Newtonian gravitational field, i.e., the Newtonian potential [/n- (The matter fields for 
a perfect fluid are the same in a "relativistic" theory as in Newtonian gravity (NG), up 
to slight modifications.) Point a) is easily fulfilled, once it is recognized d] that 
there is an exact similarity transformation in NG, which is appropriate to describe the 
weak-field limit in NG itself [IGj. This immediately suggests defining the family (Sa) by 
applying the similarity transformation of NG to the initial data defining a gravitating 
system in the investigated theory [121 CHI — preferably to that initial data, of a general- 
enough nature, which precisely defines the system of interest, S [12]. Of course, to use 
the Newtonian transformation demands that point b) has been solved, which depends 
on the precise equations of the theory. 

The application of these principles to the scalar ether-theory has been done in detail 
for its first version ^H]- (See Ref. ^7], Sect. 2, for a synopsis and a few complementary 
points.) The modifications to be done for the present version are straightforward. The 
definition of the metric (^J)-® and that of the scalar field ip ()22p imply that condition i) 
is equivalent to asking that ^/''^'^^ — > as A — > 0. Therefore, we may define a dimensionless 
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weak-field parameter simply as 

A = Sup,,M^(x) (42) 

(at the initial time, say). Moreover, from the scalar field equation we see that 

V = c^ip (43) 

satisfies the wave equation with the same r.h.s. (in the Newtonian limit where a ^ p) 
as Poisson's equation of NG, and the retardation effects should become neghgible in the 
Newtonian limit. Hence, \^ is a natural equivalent of ?7n. Thus, the Newtonian limit is 
defined by the same family of initial conditions as in the first version , though with 



the new definition ()43j) of V: at the initial time, ^ 

p(^) (x) = A2p(i) (x) , p*(^) (x) = Ap*(i) (x) , (44) 

V(^) (x) = A V(^) (x) , drV^^^ (x) = X'/^dTV'^^^ (x) , (45) 

uW(x) = yAu«(x). (46) 



The system Sa is hence defined as the solution of the above initial- value problem P^. 
One expects that the solution fields admit expansions in powers of A, whose dominant 
terms have the same orders in A as the initial conditions. It is then easy to check 
that, by adopting [M\x = A[M] and [T]^ = [T]/a/A as the new units for the system Sa 
(where [M] and [T] are the starting units of mass and time), all fields become order A°, 
and the small parameter A is proportional to (indeed A = (cq/c)^, where Cq is the 
velocity of light in the starting units). Thus, the derivation of the IPN expansions and 
expanded equations is very easy. At the first PNA, one writes first-order expansions in 
this parameter for the independent fields V^, p, u: 

V = Vo+V,/c^+0{c-'), p = po+pi/c^+0{c-^), u = uo+ui/c2+0(c-^), (47) 

and one deduces expansions for the other fields. (Of course, all fields depend on the 
small parameter A oc 1/c^.) In these varying units, we have T = X^^'^Tq where Tq is 
the "true" time, i.e., that measured in fixed units. Hence cT = cqTq is proportional 
to the true time. But since the true velocities in system Sa vary like A^''^ [Eq. (HH))], 
it is T, not cT oc Tq, which remains nearly the same, as A is varied, for one orbital 
period of a given body in the Newtonian limit. Therefore, in this limit, thus for PN 
expansions, one must take = T as the time variable, in the varying units utilized for 
the expansions. This means that, in these units, the expansions are first of all valid at 
fixed values of x and T; and one can differentiate them with respect to these variables, 
because it is reasonable to expect that the expansions are uniform in x taken in the 
"near zone" occupied by the gravitating system, and in T taken in an interval where 
the system remains quasi-periodic (16ij . 

^ Since the given system is assumed to correspond to a small value Ao <C 1, the transformation goes 
first from Aq to A = 1, and then fr'om A = 1 to the arbitrary value A. This amounts to substituting 
^ = A/Ao for A, and p'^^°\ etc., for p^^\ etc. ^H]- Moreover, we consider a barotropic fluid: p* — F{p). 
Thus, the initial conditions for p and for p* are actually not independent, and one must assume 
that F'^^^p) = XF^^^\~^p) Note that the small parameter A considered in the present paper 

corresponds to e^, where e is that used in Ref. 



14 



5.2 Main expansions and expanded equations 

Inserting (07j)i into the scalar field equation ()28|) and accounting for the fact that the 
time variable is x'^ = T (in the varying units utilized), yields after powers identification: 

AV^o = -AnGao, AVi = -AnGai + O^Vq, (48) 

where a = ao + <Ji/c^ + 0(c~^) is the IPN expansion of the active mass density. Thus, 
the retardation effect disappears in the PN expansions. (However, the "propagating" 
(hyperbolic) character of the gravitational equations is maintained through the fact 
that an initial- value problem is considered.) From (j48p with appropriate boundary 
conditions {U = 0(l/r) and gradt/ = 0(l/r^) as r — > oo) jTH], it follows that U = Vq 
is the Newtonian potential associated with cto: 



U = Vo = N.P.[ao] (n.P.[t]{X,T) = G J r(x, T)dV(x)/ |X - x|^ 



(49) 



(V will denote the Euclidean volume measure on the space M), and that (imposing the 
same boundary conditions to B as those for U) 

V,=B + -^, B^N.P.[a,], (50) 

with 

W(X,T) = G j |X-x|ao(x,r)dV(x)/2. (51) 

The mass centers will be defined as barycenters of p, the density of rest-mass in the 
preferred frame and with respect to the Euclidean volume measure V It is related 
to the proper rest-mass density p* by Lorentz and gravitational contraction [21] , so that 
p = P*7,;y^/ \f^) hence from (fT^: 

p = p*^Jf3\ (52) 
Using this and the definitions ()28|) 9 and ()35|) . and since we have from 1)49^ : 

/? = e-^/^' = l-[//c2 + 0(c-^), (53) 

we get: 



Po = Po = o-Q = 6q (54) 



and 



u 



Pi=Pt+Po( Y + 3^7), (55) 



(Tl = pi + po Y - t/ + Ho^ , 
Oi=o^+Po + 4pof/ = Pi + Po Y + + ^o) + 



(56) 



(57) 
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The expansion of the equation of motion (jHUj) and the energy equation (|H7|) gives, 
at the order zero: 

dripoUo) + djipoUQul) = PqU^ - Po,h (58) 

dTPo + d,{poui) = 0, (59) 

which are just the Newtonian equations. The expanded equations of the order one in 
1/c^ are: 

dT{pou\ + 0iui) + djipoui^ui + pou\ul + ^lU^n^) - poiv^QdrU + uldiU) = 

= aiU^i + po^i.i + PoU^i - 2f/po,i - Pi,i , (60) 



dT{wo + pi) + dj[{wo+po + pi)ul + pou{] = -podrll, Wq = PoCy + '^o - U) . (61) 

Combining , the continuity equation , and the 0-order expansion of the isentropy 
equation: 

dno = -pod(l/po), (62) 
one gets in a standard way the Newtonian energy equation: 

OtWo + dj [{wo + Po)ul] = -po OtU. (63) 

Subtracting (j63|l from (j6T| gives us 

drpi + dj{piul + pqu{) = 0, (64) 

which (together with (j59p ) means that mass is conserved at the first PNA of the scalar 
theory, also in this second version. 

5.3 Application: gravitational effects on light rays 

The effects of a gravitational field on an electromagnetic ray, seen as a "photon" (a test 
particle with zero rest-mass), represent the most practically-important modification to 
NG. In this theory, they can be obtained by applying the extension © of Newton's 
second law, in which Fq = and the mass content of the energy E = hv has to be 
substituted for the inertial mass m{v). In the new version of the scalar theory, things 
go in close parallel with the former version, based on an anisotropic space contraction 
[SZlIll: 

• i) The main step is the recognition ,37^ that the PN equation of motion for a 
photon, obtained thus, coincides with the PN expansion of the geodesic equation 
for a light-like particle in the space-time metric 7, because the Fg^ = ^g^'^dogkj 
Christoffel symbols of 7 are 0(c~^) and the T'jl = ^j^^dogjk are 0(c~^) (with 
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= T as the time coordinate). This holds true in the present version based 
on Eq. (jni) for the space metric, because the same expansion [Eq. above] 
apphes as in the former version. Therefore, to compute the effects of a weak 
gravitational field on light rays, one has to study the PN expansion of 7 in the 
relevant reference frame: that frame Ev which moves with the velocity V in the 
preferred frame, assumed constant and small as compared with c, of the mass- 
center of the gravitating system. In coordinates {x'^) that are Galilean for the 
flat metric 7° and adapted to the frame Ev, the components 7^^ of 7 are deduced 
from its components in the preferred frame (Eqs. Q and (jH))) by a special Lorentz 
transform, relative to 7° hence the 7^^, 's depend only on the field j3 (not 
any more on its derivatives, as was the case with the former version), and on the 
velocity V. 

• ii) Inserting the expansion (jK!^ of /3, one gets the PN expansion of the 'y'^^, 's. 
The PN expansion of 799 is enough to compute the gravitational redshift. It is 
still 7qq = 1 — 2Uc~^ + 0{c~^) with U the Newtonian potential: this holds true 
in the present version (in particular, Eq. (52) of Ref. [37] holds true). To get 
the other two effects of gravitation on light rays, namely the deflection and the 
time delay, one needs to compute the PN expansion of all components 7|jj,. One 
finds easily that, as before [37], 7oj = 0(c~^) (in fact 79^ = for z = 2 and 
3, now); such 7Qj = 0(c~^) component (s) have (has) no influence on the PN 
equation of motion of a photon (see the equation after Eq. (9.2.4) in Weinberg 
PS], and see Eqs. (9.1.16), (9.1.19) and (9.1.21) there). And one finds that 
7^ = —(1 + 2Uc~'^)6ij + 0(c~^). Thus, in the new version of the scalar ether- 
theory, the PN equation of motion of a photon coincides with the PN geodesic 
equation of motion of a photon in the so-called ^3] "standard PN metric" of GR, 
and this is true also in the relevant frame Ev. In particular, in the SSS case, the 
formulas for the PN effects on photons are the same as those derived from the 
(space-) isotropic form of Schwarzschild's metric — or also from the harmonic form 
of the Schwarzschild metric (which is the SSS solution of the RTG |38j), since its 
PN approximation is space-isotropic [201 coincides with the PN approximation 
of the isotropic form. These predictions are accurately confirmed by observations 



6 PN equations of motion of the mass centers 

As already mentioned, the mass centers are defined Jj'j as local barycenters of the den- 
sity of rest-mass in the preferred frame, p or rather Pexact? Eq- (|S2I)- (Henceforth, the 
index will be omitted for the zero-order (Newtonian) quantities, for conciseness; there- 
fore, the exact quantities, when needed, are denoted by the index "exact.") Integrating 
Eq. (jUUj) in the (time-dependent) domain occupied by body (a) (a = 1, A^) in the 
preferred frame E gives 




(65) 
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with 

fi = (f^i + P)U^ + pVi,^ - 2Up,^ + piu'drU + n^U,i). (66) 
Accounting for Eq. (j37j) and for Eq. (3.21) of Ref. P7j, we get: 

Ml -d\ + = + K"', (67) 
which is Eq. (4.9) of Ref. ^Zj, and with, as there, 

Ml= ! pidV, Mlsi,= [ pixdV(x), (68) 

but with modified definitions of J"* and fT"*: 

= [ [p + p(uV2 + n + 3[/)]mMV, (69) 

JDa 

/ {a,U, + pV^,)dV, (70) 



and 



^ai ^ /■ ^_2Up^^+pU^,+pu'dTU+pu^U,i]d\/ = [ [3pU^i+pu'dTU+pu'^U,i](N = Kf+K^'+K^\ 



Da 

(71) 



Together with the Newtonian equation, Eq. (j67|) allows in principle to calculate the 
IPN motion of the mass centers: due to Eq. (3.15) of Ref. ^7], the IPN acceleration 
of the mass center of body (a) is given by 

... . Mi(ai-a) , , 

in which and a are the Newtonian mass and acceleration. Equation ()67|) may be 
made tractable for celestial mechanics, as was done in Ref. JHl for the former version of 
the theory, by taking benefit of: a) the good separation between bodies, and b) the fact 
that the main celestial bodies are nearly spherical. This is left to a future work. Here, 
we will study the point-particle limit of this equation and will show that the deadly 
violation of the WEP for a small body, which was found with the former version of the 
theory does not exist any more with the new version. 



7 Point-particle limit and the WEP 

In order to define that limit rigorously and generally, we consider (as in Ref. J3]) a 
family of IPN systems, that are identical up to the size of the body numbered (1): 
this size is a small parameter ^. We have to expand as ^ ^ the integrals (j69|) . 
(frOj) . and ()71|) . for the small body, i.e. a = 1. (As to the zero-order (Newtonian) 
acceleration of the small body, it tends towards the acceleration of a test particle in 
the Newtonian field U^""^ of the other bodies 15J, as expected.) To do that, we use the 
simplifying assumption according to which the Newtonian motion of the small body is 
a rigid motion. The calculations are very similar to those ^H] with the previous version, 
though simpler for the K"-^ integral; hence, we shall be concise. 
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7.1 The general case 



The modification of tlie calculations in Ref. ^3], Sect. 3, is immediate for 1°'^ and K"-^. 
We get [reserving henceforth the letter a for the first, small body (for which a = 1 in 
fact) and using the letter b for the other, massive bodies]: 

r = {r') = Mj|a2 + 3f/('^)(a)]a + O(e^), (73) 



K1 = o{e 



K-GM.a$:^M^-;^ + 0(e 
'a — b 



= M.a'Vt/^^Ha) + 0(e^). 



(74) 
(75) 
(76) 



As to the integral J°'\ it has the same expression as before ^7], but the PN correction 
0"! to the active mass density is now given by Eq. (|^. Therefore, the expansion (3.27) 
of Ref. ^3] remains valid for but we have now: 

aa= [ cTidV = Ml + Ma [|a2 - t/('^)(a)]|T=o + O(e'), (77) 

-'Da 

Ml = M, [ia^ + 3f/('^)(a)]|r=o + O(e^), (78) 
Pa, = [ a^i^)ix^ - a^)dV(x) = Ml{a{ - a^) + O(e'). (79) 

-'Da 

Beside and Paj, Eq- (3.27) of Ref. [12] involves also all those multipoles of the 
densities p and ai that correspond to the other bodies. 

Since all of these equations contain the (Newtonian) mass Ma as a common factor, 
it follows that the IPN acceleration A" of the (mass center of the) small body (a), 
Eqs. fl67|) and (f?^ . does not depend on its mass M^. It then also follows from these 
equations (including Eq. (3.27) of Ref. [15,) that, neglecting O(^) terms in A", it 
depends only: 

• on the current Newtonian positions and velocities of all bodies: a, b, a, b; 

• on all current IPN positions: a(i),b(i); 

• on the Newtonian masses M^ of the other bodies and on their Newtonian potential 
f/^'^) (the "external" potential for (a)); 

• and still, through the external multipoles of p and ai, on the structure of the 
other bodies. 

Thus, the IPN acceleration of a freely-falling small body is independent of its mass, 
structure, and composition, in other words the WEP is satisfied at the IPN approxi- 
mation with the new version of the scalar theory, in contrast with what happened with 
the former version |15j . 
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7.2 Comparison with a test particle in the case with one SSS 
massive body 

Since the WEP is satisfied, it seems obvious that the IPN acceleration of a small body 
should be equal, at the point-particle limit ^0), to that of a test particle — and this 
in the general case. We check this in the particular case where, beside the small body 
(1), there is just one massive body (2), whose mass center stays fixed at the origin in 
the preferred frame, and whose Newtonian density p is spherically symmetric. ^ We 
note 



Ml, M = M. 



2, 



V = X, r 



n 



x/r, xi 



c^(a(i) 



Adapting Sect. 4 of Ref. [T3], we find without difficulty: 



GM 



— m- 



v2 GM 

h 3 

2 r 



8M 



(XidV = 



D2 



n + 4(v.n)v 
pf/dV 



1 

'^2 



D2 



(80) 

(81) 
(82) 



GM 



m- 



GM GM 17 e 

2 4 

r To 



3 M 
GM 



1 / / 

n H — (3(xi.n)n — Xi 



= -m- 



(ro = r\T=o), 
(83) 

(84) 



Inserting these values into Eq. (jU7j) and putting the result into (f72j) . we get the equation 
for the IPN correction Xi to the position of the mass center of the small body: 



GM 



Xl 



-v^ + 4- 



GM 



— - + 3 n + 4 v.n V 

3 M r ^ ^ r 



i5) 



On the other hand, as recalled in Subsect. 12.21 the equation of motion of a test particle 
in the scalar theory coincides, in the present static case, with the geodesic equation 
in the relevant metric — thus, here, the metric (jH)-®, specialized to the SSS case, for 
which we get from (j2Hl) and (j^ : 

/3(T,X) = e-^, (i?=|X|>r2), (86) 
with r2 the radius of the spherical body (2), and with 

M'=[ aexactdV = M + ^ + 0(c-^), (87) 

^Strictly speaking, body (2) is gravitationally influenced by the small body, hence it cannot stay 
exactly at rest in the preferred frame. However, the PN acceleration of the massive body (2), due to 
the small body (1), is 0{^^) (Ref. ^^li Sect. 2), so that we may forget this influence for the present 
purpose. 
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thus in Cartesian coordinates (X*) for the Euchdean metric g^: 

2 



GM' fGM' 
1-2^^ + 2 



cm 



cm 



0{c 



{dx 



0\2 



1 + 2^ + 0(0 



bi^dX'dX\ 



which is the SSS form of the standard PN metric of GR. The corresponding (complete) 
equation of motion is given by Weinberg |29| (Eq. (9.5.3) with here e = 0, C = 0, and 
= — GM'/r(i)). In our notation, this is 



GM' 



(1) 



-n(i) 



-^(1) + 4- 



.GM' 



n(i) +4(v(i).n(i))v(i) 



with 



X(i) = X + xi/c^ 



^9) 



(90) 



r(i) = |x(i)| , n(i) = X(i)/r(i), V(i) = X(i). 

Now, writing r(i), n(i) and V(i) as first-order expansions in c~^, and then inserting (jUUj) 
and (jS7j) into (|H^. one finds easily that the latter decomposes into an equation for the 
order 0, which is the Newtonian equation x = — ^^n, and an equation for the order 1 
in c~^, which is exactly Eq. (|85|). This proves that indeed, the IPN acceleration of a 
small body is equal, at the point-particle limit, to that of a test particle — at least in 
the SSS case. 



8 Conclusion 

The investigated theory starts from a heuristic interpretation [221 gravity as the 
pressure force exerted on the elementary particles by an universal fluid or "constitu- 
tive ether" j2Hj, of which these particles themselves would be just local organizations. 
This leads naturally to assuming that gravity affects the physical standards of space 
and time, by an analogy with the effects of a uniform motion that are at the basis of 
Lorentz-Poincare (special) relativity j22I- However, the contraction of physical objects 
in a gravitational field, as it appears in terms of the "unaffected" Euclidean metric, may 
either occur in one direction only j21], as for the Lorentz contraction, or else 1^3031221 
it may affect all (infinitesimal) directions equally. The first version of the theory was 
based on a unidirectional contraction, and passed a number of tests [21 El , but it has 
been discarded by a significant violation of the weak equivalence principle (WEP), 
which has been found to occur for extended bodies at the point-particle limit [Kj. In 
the present paper, a new version of the theory, based on a locally isotropic contraction, 
has been fully constructed. It has been shown that the new version also explains the 
gravitational effects on light rays (Subsect. 15. 3|) . Being based on the same dynamics as 
the former version, and being also based on a wave equation for the scalar gravitational 
field, it should lead, as did the former version to a "quadrupole formula" similar 
to that used in GR to analyse the data of binary pulsars |i41j. Moreover, because the 
metric in the new version is similar to the "standard PN metric" of GR while the local 
equations of motion are also similar to those of GR, the celestial mechanics of that 
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theory should improve over Newtonian celestial mechanics. These two points will 
have to be checked in a future work. 

It has been proved here that the present new version of the theory solves completely 
the problem with the WEP, that occurred at the first post-Newtonian approximation 
in the former version (Sect. [7j). When that (deadly) WEP violation had been found, 
it had been argued J3] that the reason for it was the dependence of the PN spatial 
metric on the spatial derivatives of the Newtonian potential U. By switching to an 
isotropic space metric, whose PN form depends on U but not on its derivatives, we 
indeed suppressed the WEP violation in the present new version. 
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